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ABSTRACT 

We present submillimetre data for the L1689 cloud in the p Ophiuchi molecular cloud 
complex. We detect a number of starless and prestellar cores and protostellar en- 
velopes. We also detect a number of filaments for the first time in the submillimetre 
continuum that are parallel both to each other, and to filaments observed in the neigh- 
bouring L1688 cloud. These filaments are also seen in the ^'^CO observations of L1689. 
The filaments contain all of the star-formation activity in the cloud. L1689 lies next to 
the well studied L1688 cloud that contains the p Oph-A core. L1688 has a much more 
active star-formation history than L1689 despite their apparent similarity in ^^CO 
data. Hence we label L1689 as the dog that didn't bark. We endeavour to explain 
this apparent anomaly by comparing the total mass of each cloud that is currently 
in the form of dense material such as prestellar cores. We note firstly that L1688 is 
more massive than L1689, but we also find that when normalised to the total mass of 
each cloud, the L1689 cloud has a much lower percentage of mass in dense cores than 
L1688. We attribute this to the hypothesis of iLorenl l)1989f) that the star formation 
in the p Ophiuchi complex is being affected and probably dominated by the external 
influence of the nearby Upper Scorpius OB association and predominantly by a Sco. 
L1689 is further from a Sco and is therefore less active. The influence of a Sco appears 
nonetheless to have created the fllaments that we observe in L1689. 
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1 INTRODUCTION 

The earliest stages of low-mass (O.2-2M0) star forma- 
tion arc becoming reasonably well understood (see e.g. 
[And re. Ward- Thompson, & Barsonv 2000 for a review). The 
prestellar core stage llWard- Thompson et alj Il994ll is the 
phase in which a molecular cloud core has become gravi- 
tationally bound. Thereafter gravitational collapse sets in 
and a central hydrostatic protostar forms, which is known 
as a C lass protostar (Andre, Ward-Thompson, & Barsonv 
Il993l) . Once half of the final mass has accreted onto 
the central object it is kn own as a Class I proto- 
star iWilking. Lada. fc Youn3 fT989). and it subsequently 
evolves through the Class II & III young stellar object (YSO) 
phases jLada 1987). Debate continues over the details of this 
evolutionary process, and different molecular clouds appear 
to be evolving according to a different interplay of physical 
mechanisms. These include magnetic fields, turbulence and 
feedback from previous episodes of star formation. Submil- 
limetre studies of different star-forming regions are required 
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to help clarify the evolutionary process and determine the 
effect of these mechanisms. 

The L1689 molecular cloud is part of the p Ophiuchi 
molecular cloud complex, w hich is located at a distance of 
128 ± 12 pc from the sun teertout. Robicho n. fc ArenorJ 
1999;). It was first detected by Lynds (1962) in a large study 
of the dark nebulae detected in the Palomar Observatory 
Sky Survey plates. The p Ophiuchi complex was extensively 
mapped bv iLorerl lll989h using '^^CO at a resolution of 2.4 
arcmin. The resulting map of the cloud is shown in Fig. 
The data show that the molecular cloud is composed of a 
number of sub-clouds, the most massive of which is L1688. 

The majority of studies of the p Ophiuchi molecular 
cloud to date have concentrated on the L1688 sub-cloud, 
also known as the p Oph main cloud. A number of filamen- 
tary clouds (L1709, L1740, L1744, L1755 and L1765) extend 
from L1688 in a north easterly direction, and are often called 
the streamers or the cobwebs of Ophiuchus. To the south- 
east of L1688 lies L1689, which also has filamentary clouds 
(L1712 and L 1729) extending in roughly the same direction. 
iLorenl Jl98Sl) measured masses for L1688 and L1689 of 1447 
Mq and 566 Mq respectively, while the components of the 
streamers have lower masses of 100-300 M0 each. The di- 
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Figure 1. "CO map of the p Ophiuchi complex, from iLorenl il98Hl) . The contour levels give antenna temperatures r^(^'^CO) of 4, 5, 
6, 7, 8, 10, 12, 14, 18 and 20 K. The names of the clouds are marked. 
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Figure 2. A composite map of the L1689 cloud showing the scan-mapped fields (white outline and labels). The background 0.5 ^m 
image is taken from the Digital Sky Survey, obtained using the Sky View interface. The contours are "CO emission at T^{^^CO) of 4, 
5, 6, 7, 8, 10, 12, 14, 18 and 20 K iLoren 19891 . 
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Table 1. The positions of eacii of the scan-mapped fields in the L1689 molecular cloud. The atmospheric conditions during each 
observation are indicated by the 850 ^m zenith optical depth listed in column 5. The field names correspond to those shown in Fig.l2l 



Map centre 



r lelu 


KA 


Dec. 






Name 


/ onnn\ 
(^zUUU ) 


i onnn\ 
t^zUUU ) 


U 1 JJate 


T850m™ 


R31-49 


16:31:39.6 


-24:49:43 


1999-Mar-09 


0.34 


_ 






2000-Apr-09 


0.26 


R31-58 


16:31:56.7 


-24:58:52 


1999-Mar-09 


0.38 


_ 


_ 


_ 


2000-Apr-lO 


0.22 


R32-45 


16:32:11.5 


-24:45:21 


2000-Apr-09 


0.22 


R32-03 


16:32:14.1 


-25:03:10 


1999-Aug-08 


0.22 


_ 


_ 


_ 


1999-Aug-09 


0.38 


R32-07 


16:32:22.7 


-25:07:50 


1999-Apr-08 


0.22 


_ 


_ 


_ 


2000-Apr-lO 


0.18 


R32-37 


16:32:26.1 


-24:37:20 


2000-Apr-09 


0.22 


R32-28 


16:32:27.6 


-24:28:25 


1999-Mar-09 


0.30 


R32-00 


16:32:38.8 


-25:00:49 


2000-Apr-lO 


0.18 


R33-22 


16:33:01.9 


-24:22:17 


2000-Apr-12 


0.30 


R33-37b 


16:33:02.1 


-24:37:17 


1999-Aug-09 


0.26 








2000- Apr- 11 


0.22 


R33-32 


16:33:17.0 


-24:32:16 


1999-Mar-09 


0.34 








2000-Apr-lO 


0.18 








2000-Apr-12 


0.30 


R33-46 


16:33:22.3 


-24:46:16 


2000-Apr-12 


0.30 


R33-37a 


16:33:38.1 


-24:37:15 


2000- Apr- 11 


0.22 


R34-28 


16:34:01.9 


-24:28:13 


1999-Mar-09 


0.38 








1999-Aug-08 


0.26 


R34-37 


16:34:14.1 


-24:37:13 


2000- Apr- 11 


0.22 


R34-38 


16:34:50.2 


-24:38:04 


1999-Aug-08 


0.26 








1999-Aug-09 


0.30 
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to understand the similarities and differences between L1688 
and L1689. 



2 OBSERVATIONS 

The observations were carried out using the Submillimetre 
Common User Bolometer Array (SCUBA) on the James 
Clerk Maxwell Telescope (JCMT). This instrument takes 
observations at 450 and 850 /^m simultaneously through the 
use of a dichroic beam-splitter. The telescope has a resolu- 
tion of 8 arcsec at 450 and 14 arcsec at 850 /^m. 

Observations of L1689 were carried out over several 
nights between 1999 March and 2000 April using the scan- 
map observing mode as listed in Table A scan-map is 
made by scanning the array across the sky. The scan direc- 
tion is 15.5° from the axis of the array in order to achieve 
Nyquist sampling. The array is rastered across the sky to 
build up a map several arcminutes in extent. 

Time-dependent variations in the sky emission were re- 
moved by chopping the secondary mirror at 7.8 Hz. Due 
to a scan-map being larger in size than the chop throw, 
each source in the map appears as a positive and a negative 
source. In order to remove this dual-beam function, each re- 
gion is mapped six times, w ith chop throws of 30, 44 and 68 
arcsec in both RA and Dec (|Emerso3[l993). The dual-beam 
function is removed from each map in Fourier space by di- 
viding each map by the Fourier transform of the dual-beam 



rection of the streamers turns to the south at the positions 
of both L1688 and L1689, forming a horseshoe shape in the 
case of L1688 and a boomerang shape in the case of L1689 
(see Fig. nj. 

Numerous studies of the L1688 cloud have revealed a 
rich history of star formation. Prompted by the associa- 
tion of early-type stars with the p Ophiuchi dark cloud, 
iGrasdalen. Strom, fc StromI (^T^) surveyed the cloud at 
2 /im and discovered a large number of embedded young 
stars . Subsequent observati ons at a variety of wavelengths 



e.g. Wilking fc Lada 1983,: lA ndre. Montmerle. fc Feigelson 



IWard-Thompson et al.l Il989l: IWilking^ et al 



1987 



19891: iGreene et all Il994l: ICasairova et all |199 



Motte. Andre, fc NeriT 119981: Ijohnstone et alJ 1200 



Johnstone. Di Francesco fc Kirk| |2004|) have revealed 



large numbers of objects at every stage of star formation. 
Indeed, studies of this region have greatly increased our 
knowledge of the star-formation process as a whole. 

In contrast, the L1689 cloud appears to have rela- 
tively little star-forming activity, concentrated in a small 
number of isolated sources such as IRAS 16293-2422 
jLoren. Wootten. fc Wilkind Il99d) . As such, we have la- 
belled L1689 'the dog that didn't bark'. In this paper, we 
endeavour to find out why there is such an apparent dis- 
crepancy between the star-formation activity in L1688 and 
L1689 when they appear so similar in CO maps - see Fig.0 
We report on submillimetre observations of L1689 and com- 
pare them with those taken at other wavelengths, in order 
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Figure 3. The L1689 cloud at 850 fim. The contour levels are at 3(t, 5cr, lOo", 20(t, 40a" and 80(t, where a is the mean noise level over 
the whole of the map, equal to 20 mjy/beam. Contours are based on the map smoothed to a resolution of 18 arcsec. The map regions 
in Table |5] are shown as dotted lines. 



Table 2. Table showing the noise levels for each region of the L1689 maps. The position of each of the regions is illustrated in Fig. 1^ 





Region 


Centre 


Icr noise 


(Jy/beam) 


Region Name 


RA (2000) 


Dec. (2000) 


850 Atm 


450 fim 


L 1689- West 


16:31:40 


-24:50:00 


0.029 


2.4 


L1689-South 


16:32:10 


-25:03:00 


0.016 


0.51 


L1689-NorthWest 


16:32:30 


-24:35:00 


0.047 


3.4 


L1689-Centre 


16:33:20 


-24:47:00 


0.034 


1.3 


L1689-North 


16:33:30 


-24:33:00 


0.023 


0.87 


L1689-East 


16:34:50 


-24:38:00 


0.015 


0.39 
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Figure 4. Composite map showing the 850 lira scan-map of 
L1689 as a grey-scale. The superposed contours are the 850-^m 
data smoothed to a resolution of 1 arcmin to highlight the large- 
scale structure. The names of the filaments are also marked (see 
text for details). 



function, which is a sinusoid. The multiple chop-throws al- 
low for cleaner removal of the dual beam function in Fourier 
space. The maps are them combined, weighting each map 
to minimise the noise introduced at the spatial frequencies 
that correspond to zeroes in the sinusoids. Finally the map 
is converted back into normal space, where it no longer con- 
tains the negative sources dJenncs s & Lightfoot 200^). 

The submillimetre zenith opacity at 450 and 850 /xm 
was determined using the 'skydip' method and by com- 
parison with polynomial fits to the 1.3 mm sky opacity 
data, me asured at the Caltech Submillimeter Observatory 
jArchibald et al. 2002). The average zenith optical depth at 
850 /xm is listed in Table for each observation. The tele- 
scope pointing was checked at regular intervals throughout 
the nights using planets, secondary calibrators and standard 
pointing sources. 

The data were reduced in the normal way using 
the S CUBA User Reduction Facility jjenness fc Lightfood 
l2000t) . Calibration was performed using observations of the 
planet Uranus taken during each shift. We estimate that 
the absolute calibration uncertainty is ±10% at 850 ^m and 
±30% at 450/im, based on the consistency and reproducibil- 
ity of the calibration from map to map. 



3 RESULTS 

The map of L1689 is made up of a number of scan-maps, 
each of which is 12 x 12 arcmin in size. The positions of the 
sc an-maps were chosen to trace the ^^CO emission mapped 
bv lLorenI lll989l) . This is shown graphically in Fig. |21 where 
the mapped area is overlaid in white on the black ^"^CO 




RA (2000) 

Figure 5. An enlarged view of Filament 1 at 850 fira. The sources 
SMM 1 — 16 that are associated with this filament are labelled. 

contours jLorenlll989l) . The names and approximate posi- 
tions of the different scan-map fields are shown in white. 
The backgr ound 0.5 /xm image is taken from the Digitized 
Sky Survev llLaskeJl994ri obtai ned using the Sky View inter- 
face llMcGlvnn fc ScoUicklll994^ . The mapping strategy was 
to offset each scan-map from the neighbouring maps by ap- 
proximately half a map-length. This was to avoid potential 
'striping' caused by the increased noise levels at the scan- 
map edges. The name of each scan-map is taken from the 
J2000 coordinates of the scan-map centre. The observation 
details are summarised in Tabled 

Fig. El shows the data for the L1689 cloud at 850 /im. 
The area covered is approximately 0.5 deg^, which at a dis- 
tance of 128 pc is equal to ~ 2.5 pc^. As a result of the 
different weather conditions at the telescope when the maps 
were made, and also the varying integration time per point, 
the noise level varies across the map. We have therefore split 
the map up into a number of regions, and measured the mean 
1(T noise level for each region. These are given in Table |21 
The extent of the different regions, and our names for them, 
are illustrated with dashed lines in Fig. |5] 

A number of sources are seen in the 850 fira map. We 
have labelled these SMM 1—21 in order of increasing declina- 
tion. The sources were identified by eye by selecting regions 
of flux density greater then 4-(t above the background level. 
All the sources are smaller than ~ 1 — 2 arcmin. Any struc- 
ture larger than this scale is deemed to be part of the fila- 
mentary structure of the cloud, and is discussed below. The 
identities of the relatively compact (i.e. < 1 arcmin) sources 
are discussed in the following section. The only source de- 
tected at 450 ^m is the well-st udied binary /multiple Class 
protostar IRAS 16293-2422 JWoottenlll989l: iMundv et all 
119921: IWalker. Carlstrom. fc Bieginslll993h . 
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Table 3. The position and measured flux density of each source in the L1689 scan-maps at 850 fim and 450 fim. Column 4 indicates 
whether or not the core is extended in the 14 arcsec 850 ^m JCMT beam. The peak flux density is quoted for all sources (columns 7 and 
9), with 3<J upper limits given for undetected sources. The integrated flux density within an aperture is given for extended sources in 
columns 8 and 10. The semi-major Sz semi-minor axes of these apertures are given in columns 5 and 6. The names given to these sources 
in other surveys are given in Table HI 



Aperture 850 fim 450 ^m 



Source RA Dec. Ext. Semi-major Semi-minor Peak Int. Peak Int. 

Name (2000) (2000) (arcsec) (arcsec) (Jy/beam) (Jy) (Jy/beam) (Jy) 



SMM 1 


16:32:12.5 


-25:03:53 


Y 


32 


24 


0.15 


0.75 


< 1.5 




SMM 2 


16:32:16.1 


-25:02:44 


Y 


67 


29 


0.14 


1.6 


< 1.5 




SMM 3 


16:32:02.7 


-25:01:48 


Y 


24 


22 


0.15 


0.50 


< 1.5 




SMM 4 


16:32:04.1 


-25:00:53 


Y 


38 


22 


0.18 


1.1 


< 1.5 




SMM 5 


16:32:03.0 


-24:59:48 


Y 


37 


22 


0.15 


0.67 


< 1.5 




SMM 6 


16:31:54.6 


-24:58:22 


Y 


30 


20 


0.22 


0.94 


< 1.5 




SMM 7 


16:32:06.3 


-24:58:01 


Y 


41 


27 


0.11 


0.49 


< 1.5 




SMM 8 


16:31:57.5 


-24:57:40 


Y 


43 


25 


0.30 


1.4 


< 1.5 




SMM 9 


16:32:01.0 


-24:56:44 


N 


13 


13 


0.44 




< 1.5 




SMM 10 


16:32:08.6 


-24:56:41 


Y 


39 


27 


0.11 


0.45 


< 1.5 




SMM 11 


16:31:52.2 


-24:56:13 


Y 


17 


15 


0.29 


0.56 


< 1.5 




SMM 12 


16:31:53.7 


-24:55:54 


Y 


36 


21 


0.32 


1.4 


< 1.5 




SMM 13 


16:32:00.3 


-24:55:53 


Y 


31 


21 


0.23 


0.72 


< 1.5 




SMM 14 


16:31:37.5 


-24:51:09 


Y 


36 


26 


0.27 


1.5 


< 7.1 




SMM 15 


16:31:46.8 


-24:51:10 


Y 


34 


27 


0.20 


1.2 


< 7.1 




SMM 16 


16:31:39.2 


-24:49:48 


Y 


63 


51 


0.40 


7.8 


< 7.1 




SMM 17 


16:33:55.7 


-24:42:06 


N 


12 


11 


0.21 




< 2.6 




SMM 18 


16:34:48.4 


-24:38:04 


Y 


43 


30 


0.24 


2.0 


< 1.2 




SMM 19 


16:32:29.3 


-24:29:13 


Y 


48 


33 


1.64 


12 


< 10 




SMM 20 


16:32:23.0 


-24:28:40 


Y 


39 


28 


15.9 


28 


80 


180 


SMM 21 


16:33:06.2 


-24:28:38 


Y 


27 


19 


0.12 


0.36 


< 2.6 





Table |3| gives the measured flux density for each source 
detected in the map. The peak flux density is given in 
Jy/beam for each source at both 450 and 850 fim. 3cr upper 
limits are given for sources that are undetected at 450 fim. 
The integrated flux density is also given for sources that are 
extended in the 14 arcsec JCMT beam. The semi-major and 
semi-minor axes of the elliptical apertures used to measure 
these integrated flux densities are given in columns 5 and 6 
of Table El 

Taking the map as a whole, we see that the cloud ap- 
pears to be composed of a number of approximately parallel 
filaments, that run the length of the map. These are high- 
lighted in Fig. |11 where the contours are the data after it 
has been smoothed to a resolution of 1 arcmin to illustrate 
the large scale structure. A number of sources associated 
with these filaments are also detected. The nature of these 
sources is discussed below. 

Filament 1 lies at the southern end of the map, and 
extends northwest-southeast (NW-SE) through the L1689- 
West and L1689-South regions. The filament appears to have 
fragmented into a number of objects, forming more than 
one object across its width in some cases. The degree of 
fragmentation appears to be greater at the north of the fil- 
ament, where the majority of the sources detected in this 
filament are found. The sources at the southern end of the 
filament are also seen to be less centrally condensed. Fig. |K| 
shows this filament in more detail. The 16 sources in the 
filament have been labelled. Most of these objects are sig- 
nificantly extended in the JCMT beam, and have a fairly 



complex structure. This filament contains the majority of 
the sources that are seen in the L1689 molecular cloud. 

Filament 2 runs roughly parallel to Filament 1 but 
slightly to the north, also in a NW-SE direction. There are 
no detected sources associated with this filament. The fila- 
ment has been independently detected on a number of scan- 
maps. Unfortunately, Filament 2 has not been covered very 
well by the scan-mapping strategy. There is an area that has 
not been mapped in the middle of the filament, and most 
of the area that has been mapped has been covered only 
once. In addition, the filament has mostly been caught on 
the edges of the scan-maps, where the noise levels can be of 
order 1.5x higher. As a result of these factors, the level of 
noise along the filament is higher than average. We there- 
fore suggest that this filament be re-mapped to confirm our 
findings. 

Filament 3 extends from the L1689-NorthWest region, 
in a NW-SE direction also roughly parallel to Filament 1. 
The brightest source in this filament is the IRAS 16293 bi- 
nary/multiple Class protostar (labelled here SMM 20). 
There is also a weaker source 90 arcsec (approximately 
14000 AU) to the southeast of IRAS 16293, named 16293E 
(SMM 19). It is interesting to note that the extended emis- 
sion surrounding IRAS 16293 lies parallel to the axes of the 
other filaments in the map. The filament is only detected 
in the northwest of the image, though if the line of the fll- 
ament is extrapolated, a weak unresolved source (SMM 17) 
is detected at RA (2000) = 16''33""55.7% Dec. (2000) = 
-24°42'06". 
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Figure 6. Composite map sliowing the 850 /im scan-map of L1689 as a grey-scale. Tlie dark contours are 850-/jm data smootlied to a 
resolution of 1 arcmin to highlight the large-scale structure. The light contours are ^^CO emission with Vslr of (a) 4.68 — 5.02 kms~^ 
and (b) 3.32 — 3.66 kms~^ ILoreniil989.) . and give antenna temperatures T^(^^CO) of 2, 4, 6, 8, 10 and 12 K. The continuum filaments 
can be seen at different velocities in the '^^CO data. 



Filament 4 extends from the L1689-East region at 
the east of the map, in a NW-SE direction. Like Fila- 
ment 2, it skirts along the edge of the map, therefore 
the noise is higher than average. The only object de- 
tecte d in the filament is the L1689B jMy ers fc BensonI 
^83)^prcstcllar core llWard-Thompson et al.l 11994 
Kirk. Wa rd-Tho mpson, fc Andrei |2005^ at RA (2000) 
= 16''34'"48.6% Dec. (2000) = -24°38'00", that we label 
here SMM 18. Filament 4 is the least completely mapped 
of the four filaments, but it is seen in more than one 
overlapping scan-map. We also note that the envelope of 
L1689B is elongated parallel with the filament axis. 

When observed with molecular line tracer s, the p Qph i- 
uchi cloud is also filamentary in appearance jLorenlll989l) . 
This filamentary structure is not restricted to the streamers 
that extend to the northeast, but is also s een in the larger 
L168 8 and L1689 sub-clouds JWilking fc Lada 1983: L oreni 

Iliii). 

To illustrate this, Figs|SJa)fc(b) show our 850 fim data 
smoothed to highlight the fila ments as in Fig.El ^^CO maps 
are overlaid as hght contours jLorenl ll989'). Figs. |ni;a)fc(b) 
show the ^^CO radial velocity intervals 4.68 — 5.02 kms~^ 
and 3.32 — 3.66 kms~^ respectively. As can be seen. Fila- 
ment 1 is clearly detected in ^^CO at both of these velocity 
intervals. Filaments 3 and 4 and the northern half of Fila- 
ment 2 are seen at velocities 3.32 - 3.66 kms~^ (Fig-int)- 



4 INDIVIDUAL SOURCES 

4.1 Identities of the detected sources 

In identifying the nature of submillimetre detected sources, 
there are a number of factors to take into account. To differ- 



entiate between a prestellar core and a Class protostar, one 
can look for a centimetre radio source or a coUimated molec- 
ular outflow, which would indicate the presence of a cen- 
tral hydrostatic object. Class protostars are distinguished 
from more evolved classes of protostar by a large ratio of 
submillimetre to bolometric luminosity, or the presence of a 
centrally peaked and extended envelope |,An drc et al. 200(J). 
More evolved protostars (Class I and later) can be seen in 
the infrared. In this section, we discuss the identities of each 
of the sources, and summarise the source parameters in Ta- 
ble H 

In order to estimate whether or not the sources are grav- 
itationally bound, we compare their m asses to the Bonnor- 
Ebert critical mass Mbe = 2ARal/G jBonnoJll95ilEberti 
^55), where R is the radius of the source (taken to be the ge- 
ometric mean of the semi-major and semi-minor axes given 
in Table |3J, as is the isothermal sound speed, and G is the 
gravitational constant. 

The source SMM 1 corresponds to the prestellar core 
L1689A, w hich has previou sly been mapped in the sub- 
millim ctrc (Kirk et al. 2005 ^ and also in the far-infrared 
JWard- Thompson. Andre, fc Kirkll2002^ . This core's posi- 
tion, close to the leading edge of the L1689 cloud nearest the 
Upper Scorpius OB association (see Section [5.31 . could po- 
tentially explain the unusual temperatu re pr ofile m easured 
in the far- infrared bv lWard- Thompson et al] I2OO2I) . 

The sources SMM 2-5 are all fairly extended cores that 
are not very centrally peaked. Their masses are comparable 
with the Bonnor-Ebert critical mass, as given in Table |31 
These sources have not been observed previously, with the 
exce ption of SMM 5 which w as detected in the millime- 
tre llAndre fc Montmerld Il993) . and possibly corresponds 
to the source VLA 16289, though the nature of this ob- 
ject is unknown. SMM 7 is also extended, though its mass 
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Table 4. The prestellar core masses and protostellar envelope masses for the SCUBA detections in L1689. The assumed values of the 
temperature and the dust opacity for each of the sources is listed in columns 3 & 4 (see text for details). The Bonnor-Ebert critical 
mass for the starless cores is given in column 6. The region of the cloud that each source is associated with is incl uded in column 7 . 
*IRAS 16293 is known to be a binary /mul ti ple system. References in co l umn 8 are as follows : ^Bacmann ct al. 2000, ^ Green^e^aL^^94, 
- Bontemps et alJIgOOll iKirk et alTl2005l ^Andre fc Montmerlj 11994 ' ^Mundv et"ai]|l992l '^Castets et al. 2001, ^Martin et alj|l998l 
^Andre et al.lll996L 



Source 


Class 


Temp 




Mass 


Mbe 


"CO 


Other 


Name 




(K) 


(cm^g^-"-) 


(Mq) 


(Mq) 


region 


Names 


SMM 1 


Prestellar 


12 


0.01 


0.3 


0.4 


R55 


L1689A'' 


SMM 2 


Prestellar 


12 


0.01 


0.6 


0.7 


R55 




SMM 3 


Prestellar 


12 


0.01 


0.2 


0.3 


R55 




SMM 4 


Prestellar 


12 


0.01 


0.4 


0.4 


R55 




SMM 5 


Prestellar 


12 


0.01 


0.3 


0.4 


R55 




SMM 6 


Prestellar 


12 


0.01 


0.4 


0.4 


R55 




SMM 7 


Starless 


12 


0.01 


0.2 


0.5 


R55 




SMM 8 


Prestellar 


12 


0.01 


0.6 


0.5 


R55 


L1689SMM4 


SMM 9 


I 


30 


0.02 


0.02 




R55 


IRAS 16289-2450, L1689-IRS62, 
ISO-2093, L1689S-IRS674 


SMM 10 


Starless 


12 


0.01 


0.2 


0.5 


R55 




SMM 11 


II 


30 


0.03 


0.02 




R55 


L1689-IRS52, ISO-2043 


SMM 12 


Prestellar 


12 


0.01 


0.6 


0.4 


R55 




SMM 13 


Prestellar 


12 


0.01 


0.3 


0.4 


R55 




SMM 14 


Prestellar 


12 


0.01 


0.6 


0.5 


R53 




SMM 15 


Prestellar 


12 


0.01 


0.5 


0.5 


R53 


R53I 


SMM 16 


Prestellar 


12 


0.01 


3.0 


0.8 


R53 


R53I 


SMM 17 


II 


30 


0.03 


0.01 




R65 


RX J1633.9-2442** 


SMM 18 


Prestellar 


12 


0.01 


0.8 


0.5 


R65 


L1689B'' 


SMM 19 





30 


0.02 


0.6 




R57 


IRAS 16293-2422E^ 


SMM 20 


0* 


30 


0.02 


1.3 




R57 


IRAS 16293-2422'5 


SMM 21 


Starless 


12 


0.01 


0.1 


0.3 


R57 





is less than half of its Bonnor-Ebert critical mass. Therefore 
this source may well not be bound, but a transient object. 
Hence we label it as starless rather than prestellar in Ta- 
ble 4. This follows the distin ction of definitions made by 
IWard- Thompson et alJ ll200d) . who stated that the differ- 
ence between a starless core and a prestellar core is that a 
prestellar core should be consistent with being gravitation- 
ally bound. 

SMM 8 is the pre s tellar core L1689SMM, which was 
mapped by iKirk et all (|200^ using SCUBA in jiggle-map 
mode. There is no evidence of a protostar at the centre. 
The core is elongated in a north-south direction, though 
shows no evidence of fragmenting along its length. Both the 
morphology and flux density of our data are consistent with 
the earlier jiggle-map. There is a weaker source (SMM 6) to 
the south-west of SMM 8, which partially overlaps it. 

SMM 9 is the Class I protostar IRAS 16289-2450. It 
has also been detected at a n umber of wavelength s from the 
near-infrared (L1689-IRS6 -iGreene et al."l994!). the mid- 
infrared (ISO-209 - Bontein^t al.ll20 01) and the millime- 
tre (L1689S-IRS67- lAndre fc Montmerldil994ll) . This object 
is centrally condensed in the submillimetre and no extended 
emission is detected. 

There is an abrupt gap in the filament between SMM 8 
and the group of objects to the north (SMM 11, 12 and 13). 
This could be caused by the gravitational collapse of the 
filament onto the surrounding cores, or the cavity could have 
been cleared by outflows from IRAS 16289-2450 (SMM 9), 
which lies within the cavity. 



SMMIO is similar in appearance to its neighbour, 
SMM 7, and also appears not to be bound. SMM 11 is a 
Class II protostar which has p reviously been ident ified in the 
near-infrared (L168 9-IRS5 - iGreene et al.lll994l) and mid- 
infrared (ISO-204 - iBontemps et alJl20oJl . Approximately 
25 arcsec north of this (0.015 pc) lies SMM 12, which has 
not been previously observed at other wavelengths. This ob- 
ject is assumed to be prestellar in nature. SMM 13 has also 
not been detected at other wavelengths. Its mass is ~70% 
of the critical Bonnor-Ebert mass, therefore it may not be 
bound. Again we label this source as starless in Table 4. 

SMM 16 lies at the northern end of Filament 1, and is 
one of the brighter sources in the map. It appears to have 
some internal structure, though this is poorly resolved in 
the 14 arcsec beam. This suggests that SMM 16 could be 
composed of multiple unresolved sources, surrounded by a 
common envelope. The diameter of the source is approxi- 
mately 2 arcmin (0.08 pc at a distance o f 130 pc). A VLA 
survey has been conducted of this region llStine et alJll98d) . 
However, no centimetre sources coincident with SMM 16 
were detected. This indicates that the source (or sources) 
may be prestellar in nature. 

Two sources (SMM 14 & 15) are located 1-2 arcmin 
from the centre of SMM 16. Again there is no evidence of as- 
sociated protostars liStine et a l. 1988), we therefore assume 
that these objects are prestellar, based on their Bonnor- 
Ebert masses. The three sources SMM 14, 15 and 16 are co- 
incident with the starless core R53 which has been obser ved 
in absorption in the mid-infrared (iBacmann et alj|2000l) . 
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SMM 17 is a very weak isolated unresolved source, 
which lies 18 arcsec from the Class II YSO RX J1633.9-2442, 
which was identified through optical spectroscopy foUow- 
up observations to a ROSAT X-ray survey of p Ophiuchi 
iMartin et al.lll998l) . It is assumed that the submillimetre 
emission is associated with this source. 

SMM 18 is the well studied prestellar core L1689B (e.g. 
iKirk et a l."2005'. and references therein) . This was one of the 
original starless ammonia cores mapped bv iMvers fc BensonI 
il983r and subsequently determine d to be prestellar in 
nature iWard-Thompson et al.lfl994) . Subsequent mapping 
with the 1.3 mm MPIfR bolometer array at the IRAM 
30m telescope con firmed the results of the earlier survey 
i Andre et al.lll99^ , and provided a more detailed picture of 
the density profile of the core. The morphology of L1689B, 
as shown in Fig.|3 is elongated in a roughly east- west direc- 
tion. This is consistent with the 1.3 mm map ( Andre ot al. 
1996) and als o a SCUBA jiggle-map of the core l|Kirk et alH 
2005' ). Ward-Thompson et al.l i2002l) performed a greybody 
fit to the SED of L1689B between 90 ^m and 1.3 mm using 
ISO, JCMT and IRAM data, and found a best-fit tempera- 
ture of ll.^tll K. 

The brightest source in the map is SMM 20, which 
is the binary /multiple Class protostar IRAS 16293. 

This object has been well stud ie d at a num ber of 

(e.g. 'Wootton' 'l989'; 'Mundv et al.l Il992 : 



M 



«850-B850,T ' 



(1) 



i]l22^ 

200A . 



wavelen gths 

[ Walker. Carlstrom. fc B iegin g ,,1S)9! 3,; ^ Ccccarclli et al 
Castets et al ri200 O Griffin, fc Saracenol \ 

and is used as a secondary cali brator f or su bmillime- 
tre observations using SCUBA (ISandelH It ap- 

pears to be a hierarchical system, with two components 
(16293A1 & A2) separated by 47 AU, and a third source 
(1629 3B) located ap proximately 800 AU to the north- 
west llWoottenlll989h . The objects are Class protostars 
and are the sources of a quadrupolar molecular outfiow, 
though 1629 3B shows no evidence of cur rent outfiow ac- 
tivity (Walker. Carlstrom. fc Bieginelll993) . 16293A and B 
are separated by 5 arcsec and are therefore unresolved 
in the 14 arcsec JCMT beam. Molecular line studies of 
speci e s that adsorb to gra in surfaces (e.g. iCastets et al. 
I2OOII : IWakelam et al.l l2004^ indicate that the envelope of 
IRAS 16293 is composed of a hot core, approximately 2 arc- 
sec (260 AU) across, at a temperature of ~ 100 K. The core 
is surrounded by a cooler outer envelope at a temperature 
of 20 - 30 K. 

SMM 19 corresponds to 16293E, which is located ap- 
proximately 85 arcs ec (0.05 pc) to the east of IRAS 16293. 
ICastets et ail have argued that 16293E is a Class 

protostar, on the grounds of a detected molecular outflow. 
[star k et al. (2004) did not detect the outflow from 16293E, 
and argued that the source is prestellar. This may have been 
due to a sensitivity effect, so we treat this source as a Class 0. 

SMM 21 is a weak source with no previous detections 
at other wavelengths. Comparison with its Bonnor-Ebert 
critical mass indicates that it is probably not bound. We 
label it as a starless core. 



4.2 Source Masses 

Table m summarises the results of the previous section and 
gives mass estimates for all of the detected sources. The 
envelope mass M is calculated using: 



where Sgso is the 850 p,ra ffux density, D is the distance to 
the source, Kgso is the mass opacity of the gas and dust, and 
Bsbo.T is the Planck function at temperature T. 

For the starless and prestellar cores in this sample, 
we assume a value for Ki.smm of 0.005 cm^ g ~^, a s rec- 
ommended by iHenning. Michel, fc Stoenienkd 1^9^, and 
scale this to 850 fim using a canonical value of P of 1.5 
(jAndrc & Montmcrlc 1994), giving ussoum = 0.01 cm^g~^. 
These values of ft assume a gas-to-dust ratio of 100 
jHildebrandlll983f) . For the denser envelopes around Class 
and I protostars, it is assumed that the formation of 
icy mantles will increase the n i.zmm by a factor ~ 2 
(Andre & Montmerlc 1994; Hcn ning et alj ri995'l. Again, a 
canonical f3 of 1.5 is used, giving a value of Kgso/jm ~ 
0.02 cm^g~^ for these classes of object. For the higher den- 
sity circumstellar disks surrounding Class II YSOs, grain ag- 
giomeration will increase fti.smm still further jHenning et alJ 
.1995»') . .Beckwith et a l. (1990) found that the best fits to 
the submillimetre spectra of T-Tauri stars are given by 
ni.'imm ~ 0.02 cm^g~^ and /3 = 1. These values correspond 
to K850m™ = 0.03 cm^g^V 

The temperature of 1 2 K calculated for L1689B 
JWard- Thompson et alj2002D was used for all of the prestel- 
lar and starless cores. For protostellar envelopes, the vol- 
ume averaged dust tempe rature is assumed to be 30 K 
jAndro fc Montmeri3ll994l) . 



5 COMPARISON OF L1689 WITH L1688 
5.1 Cloud constituents 

It is interesting to compare the L1688 and L1689 clouds, 
as they share the same environment, yet previous studies 
have indicated that L1688 has a much higher level of star 
formation fe.g. iLoren et al]ll99Cf ). iLore n (1989) measured 
the mass of L1688 and L1689 from "CO maps to be 1447 
and 566 Mq respectively, differing by a factor of 2.6. The 
area on the sky of L1688 (measured at the 6 K contour of 
the ^^CO maps — see Fig.0 is 1.4 times larger than L1689. 
Assuming that the two clouds have approximately the same 
aspect ratio (which is a fair assumption, given their similar 
morphology), the volume of L1688 is 1.7 times that of L1689. 
The average density of L1688 is therefore approximately 1.5 
times larger than that of L1689. As a result of this, we might 
expect star formation in L1688 to proceed more rapidly than 
in L 1689. 

iLorenl 1^8^ divided the two clouds into 89 different 
regions, which have a clear spatial or velocity separation. 
These regions are used in the following analysis to compare 
the star-formation activity across the clouds. Fig. [7| shows 
the region boundaries overlaid on the ^^CO contour map. 
The names of the regions that have been mapped in the 
sub millimetre ar e mark ed on Fig. |7| (see below) . 

iLoren et alJ lll99Cfl investigated differences in the star- 
formation rates of the different regions of the p Ophiuchi 
cloud. To do this, they compared the total luminosity (L) 
of the young stars associated with each region, with the 
mass (M) of the gas, determined using the ^^CO mea- 
surements iLorenlll989ll . They measured values of L/M of 
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18"30 



RA(1950) 



Figure 7. The different regions of L1688 an d L1689, defined by 
spatial or velocity separations iLorenlllQSgl) . The region names 
are marked for areas of the clouds which have been mapped in 
the submillimetre. 



5.2 Lq/Mq for L1688 and 0.1 Lq/Mq for L1689, indicating 
a much larger star-formation efficiency (SFE) in L1688. This 
meth od of estimat i ng the SFE is biased towards high mass 
stars. iLoren et al.l (|l99ff l also calculated the ratio N/M, 
where A'^ is the number of young stars associated with each 
region. This estimate of the SFE is biased towards low mass 
stars. They again found that of all of the components of 
p Ophiuchi, L1688 has the highest SFE, with a value of 
N/M of 0.06. However, the value of N/M of 0.05 for L1689 
is only slightly lower. 

These methods are both probes of the past star- 
formation activity in the two clouds. In the following sec- 
tion, we look at the mass of the material currently involved 
in star formation within each region of the two clouds, in 
order to compare the current star-formation activity. 



5.2 Star- formation efficiency 

IJohnstone et al.l i2000l) mapped the densest parts of the 
L1688 cloud using the SCUBA camera on the JCMT. To- 
gether with the data presented here, these make an ideal 
dataset to c ompare the curren t star- formation activity of the 
two clouds. IJohnstone et al.l smoothed their map of L1688 
by twice the largest chop throw (130 arcsec), and subtracted 
this from the map. This was done in order to remove all large 
scale structure, some of which could result from residual at- 
mospheric noise or incorrectly fitted baselines in the data 
reduction process. 

For the purposes of directly comparing the two datasets, 
we carried out the same step, and subtracted the large scale 
structure. In the following analysis, all the core masses are 
calculated assuming a distance to the region of 130 pc, and 
a dust temperature of 20 K. Again, this is to allow a direct 
comparison of the two data-sets. 

The sensitivities of the two surveys are not significantly 
different. The SCUBA map of L1688 has a la sensiti vity of 
approximately 10 mjy/beam jjohnstone et al.]l2000l) . com- 
pared to our map of L1689, which has a la sensitivity of 15 
- 30 mJy/beam (see Table I^J. 

The SCUBA map of L1689 presented here covers most 



Table 5. The masses measured using different tracers for the 
L1688 and L1 689 clouds. T he mass of each ^"^CO region is given 
in column 3 <Lorenlll98gl) . The sum of the masses of all star- 
less and prestellar cores in each of the regions is given in col- 
umn 4 and simply labelled A/pro- These were calculated using 
the S CUBA data from this study for L1689, and SCUBA data 
from [ Johnst one et alJ l2000t) . as well as IRAM 30m data from 
iMotte e^aTT^^g^^ lfor L1688. All the masses quoted in the table 
have been scaled such that they correspond to an assumed dis- 
tance of 130 pc and temperature of 20 K. For each ^'^CO region, 
the ratio of the prestellar and starless mass to the total mass of 
the region is given in column 5. The prestellar mass for R26 is 
calculated using table 2 of Motte et al. 1 1998). Note the difference 
of a factor of two between the total Mprc/MisQQ for L1688 and 
L1689. 



Cloud 


"CO 


Ml 3 CO 


Mpro 




Name 


Region 


(Ms) 


(Mo) 


Mi3co 


L1688 


R22 


557 


16.4 


0.029 




R25 


243 


0.6 


0.002 




R26 


136 


0.2 


0.001 


Total 




936 


17.2 


0.018 


L1689 


R53 


40 


1.3 


0.033 




R55 


55 


1.4 


0.025 




R57 


133 


0.1 


0.0008 




R58 


1 










R59 


15 










R60 


20 










R61 


7 










R62 


7 










R63 


13 










R64 


4 










R65 


26 


0.3 


0.012 


Total 




321 




0.010 



of the "CO regions named by iLorenl (^S^ between R53 
and R65 (see Table |3 and Figure |7|l. The area covered 
by the SCUBA map of L1688 is approximately coincident 
with R22. The SCUBA map also c overs the portion of R25 
that is coincident with p-Oph F (see lLoren fc Woottenll9'8^ : 
iLoren et aT]ll990l). Though this region was not completely 
mapped by Johnston^^^alJ j20Cld) , a wider area survey of p 
Ophiuchus ljohngton^^^l] |20^i has indicated that there 
is very little star-formation activity in R25 except that con- 
tained in p-Oph F. 

The SCUBA map of L1688 does not include the p- 
Oph D core in R26, therefore data obtained using the IRAM 
30m telescope (iMotte et alJll99^ were used to extend the 
following analysis to this region. iMotte et alJ separated the 
IRAM data into large and small scale structure using a 
wavelet analysis. We use the core masses for their sma ll scale 
' clump s' (see their table 2), which, like the Johnsto ne et alJ 
l|2000|) data, have effectively had the larger scale structure 
removed. For consistency with the other data-sets, we also 
scale the R26 core masses to a temperature of 20 K and a 
distance of 130 pc. 

Table 1^ gives the mas s of each cl oud region, measured 
from the ^'^CO emission llLorenl ll989'l. The total mass of 
prestellar and starless cores within each region (Mpre) is 
given in column 4. This is calculated fro m the submillimetre 
dust emission from the cores in L1688 (iMotte et alJll998l : 
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I Johnstone et al.ll2000l) and L1689. The cores that contain 
protostellar sources have not been included in these values. 
This is so we can evaluate and compare the potential yield 
of the next generation of star formation in the two clouds. 

TableOshows that the current star-formation activity is 
highly variable across both clouds. R22 has the largest mass 
of dense cores, with 16 compared to < 2 M© in each of 
the regions of L1689. This may be reasonably expected, as 
the R22 region is significantly more massive than the other 
regions. In order to compare the fraction of the mass in dense 
cores in the different regions directly, we divide the mass in 
cores b y the mass of the region, measured from ^"^CO maps 
(Lorcn 198?). This is given for each region in column 5 of 
Table El 

When the two clouds (L1688 and L1689) are considered 
as a whole, L1688 has a much higher fraction (by a factor 
of 2) of its mass contained in dense cores, with a value of 
Mpre/Mi3co of 0.018, compared to 0.010 for L1689. Even if 
we ignore the regions of L1689 that contain no submillime- 
tre cores, and only consider R53, R55, R57 & R65, the total 
prestellar core mass ratio for L1689 only rises to 0.012. This 
is still only two-thirds of the value seen in L1688. This indi- 
cates that L1688 will probably have a higher star formation 
efficiency in the next generation of star formation. 

However, the values of Mpre/MisQQ for the most active 
regions of the two clouds (i.e. R22 in L1688 and R53 & 
R55, which make up Filament 1 in L1689) are in fact similar 
(although the absolute masses are much smaller in R53 and 
R55). This could be indicating that parts of L1689 may not 
be as infertile as was previously thought. We note that these 
three regions are located at the western edges of L1688 and 
L1689. We consider the significance of this in the following 
section. 

We finally note that the fracti on of prest ellar materia l 
in p-Oph A fsee lLoren fc WootteDlll98 6: Lor en et al.lll990ll 
is an order of ma gnitude higher tha n in the rest of R22 
(in agreement with lMotte et aLlll998f) . To calculate this, we 
compared t he mass of compact sub millimetre condensations 
in p-Oph A jjohns tone et aL||200C ) with the mass measured 
using C^^O fsee TWilking fcLadalll983l their table 2 and 
fig. 4). 

To summarise this section, we have found that on aver- 
age, the L1688 cloud has a much higher fraction of its mass 
in dense cores than L1689. The western edges of both L1688 
(R22) and L1689 (filament 1) have comparable fractions, p- 
Oph A is significantly more efficient at converting its mass 
into dense cores than any other component of the molecu- 
lar cloud, and it appears likely that the existing young star 
cluster will continue to grow in richness. The regions to the 
east of both clouds have a much lower mass fraction in dense 
cores. 



mation, which propagates through the cloud. This scenario 
has previous ly been used to explain the star form ation in 
p Oph (e.g. iLoren fc WoottenI Il98el : iLorenl Il989l) . In this 
case, we believe that the star formation in both L1688 and 
L1689 is being influenced by members of the Upper Scorpius 
OB association. 

The most massive and luminous nearby component of 
the Upper Scorpius OB association is a Sco. This is a hier- 
archical multiple containing an 09V star and a B2 III star 
and at least two other B stars (e.g. |Pigulsk^ll992^ . o Sco 
lies at a distance of 100 arcmin from L1688 (projected dis- 
tance 4 pc), and 150 arcmin from L1689 (projected distance 
6 pc). It appears to be a prime candidate for the triggering 
mechanism in both of these clouds. 

This hypothesis is consistent with the observation that 
the majority of the star formation in both clouds is occurring 
primarily in regions R22, R53 & R55, which make up paral- 
lel filaments at the edge of each cloud closest to a Sco. These 
filaments are perpendicular to the line of sight joining each 
cloud and a Sco. In addition, the extensions of both L1688 
and L1 689 to the north-east lie on a line directly away from 
a Sco jLorenl lT989). In this hypothesis, the eastern regions 
of the two clouds are shielded from the triggering source by 
the two filaments. We note that R57 is unusual in that it is 
located on the western edge of L1689, yet has a very small 
value of Mpre/Mi3co- However, R57 has recently formed the 
protostellar sources SMM 19 & 20, and if the mass associ- 
ated with these sources is included, then the fertility of R57 
becomes comparable with the regions of Filament 1. 

The influence of cr Sco has clearly had a different effect 
on the two clouds, as is evidenced by the young proto-cluster 
that has formed in L1688, and the lack of such a cluster in 
L1689. We hypothesise that this difi'erence in star-formation 
activity in the two clouds is simply due to the fact that 
L1688 is closer to a Sco. Hence it suffers a greater effect 
according to the square of the distance difference. 

The sequential star formation in this region can be 
traced back a number of generations. The Upper Scorpius 
OB association is adjacent to an older OB association called 
Upper Centaurus-Lupus. Prcibisch & Zinncckcr (2001) sug- 
gest that the massive stars in Upper Centaurus-Lupus trig- 
gered the star formation in Upper Scorpius, just as Upper 
Scorpius is triggering star formation in p Oph. It appears 
though that this is the end of the line for sequential star for- 
mation, as the only massive sta r to have formed in L1688 or 
L168 9 is the young B3 star SI llGrasdalen. Strom, fc StromI 
Il973h . It is therefore unlikely that this star will have the en- 
ergy required to propagate the star formation as effectively 
through the remainder of the p Ophiuchi cloud in the same 
way. 



5.3 Sequential star formation 

One possible explanation for the picture painted in the pre- 
vious section is that the star formation in both of the clouds 
is being affected by an external influence, such as the by- 
products of nearby young massive stars. This scenario is 
known as sequen tial or triggered star form ation, and was 
flrst suggested bv lElmegreen fc Ladal (^^. 

In this picture, the stars which have already formed, in- 
teract with the molecular cloud, triggering further star for- 



6 CONCLUSIONS 

We have mapped the L1689 cloud at submillimetre wave- 
lengths, and detected a number of cores and protostellar en- 
velopes, some of which are reported here for the first time. 
We have detected filaments in the submillimetre continuum 
that are also seen in the ^^CO observations of the cloud. All 
of the detected star-formation activity in L1689 is contained 
within these filaments, the majority being in Filament 1. 
We have compared the potential future rate of star for- 
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mation of L1689 with that of L1688 by looking at the frac- 
tion of mass of each cloud that is currently in starless and 
prestellar cores. We find that when normalised to the total 
mass of each cloud, the L1688 cloud has on average a much 
higher percentage of mass in cores. This is consistent with 
the higher star-formation activity observed in L1688. 

We deduce that L1688 and L1689 are both examples 
of triggered star formation caused by a Sco and the reason 
that the L1689 dog has not barked is that it is less massive 
than L1688 and is further away from a Sco. 
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